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We theoretically analyze weakly attenuated electromagnetic waves in quasi-two-dimensional 
(Q2D) metals in high magnetic fields. Within the chosen geometry, the magnetic field is directed 
perpendicularly to the conducting layers of a Q2D conductor. We showed that longitudinal collec- 
tive modes could propagate along the magnetic field provided that the Fermi surface is moderately 
corrugated. The considered waves speeds strongly depend on the magnetic field magnitude. Also, 
we analyzed interactions of these quantum waves with sound waves of fitting polarization and prop- 
agation direction, and we showed that such interaction may bring significant changes to the low 
temperature magnetoacoustic response of Q2D conductors. 
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I. Introduction 

As known for a long time, electromagnetic waves inci- 
dent at a surface of a metal from outside, cannot pene- 
trate into the bulk of material deeper than a thin layer 
adjoining the surface (skin layer) The suppression 
of the electric field inside the metal is caused by the re- 
sponse of conduction electrons. It occurs when the inci- 
dent wave frequency ui is smaller than the charge carriers 
plasma frequency loq, which determines the characteris- 
tic time of their response to the electromagnetic distur- 
bance. When uj > loq, the conduction electrons do not 
have sufficient time to respond, so the metal becomes 
transparent for electromagnetic waves of high frequen- 
cies. Here, we concentrate on the waves whose frequen- 
cies belong to the range r" 1 <w<u where r _1 is 
the characteristic frequency of the conduction electrons 
scattering by various types of scatterers (such as phonons 
and/or lattice imperfections). When the wave frequency 
remains within the considered range, some charge car- 
riers move along the direction of the wave propagation 
with the speed coinciding with the phase velocity of the 
wave. At weak scattering (lot > 1) such quasiparticles 
may stay "in phase" with the wave during several wave 
periods. All this time they see the electric field of the 
same magnitude and direction, and the field accelerates 
them. As a result, the electric field energy is transferred 
to the conduction electrons thus weakening the field. 

An external magnetic field B applied to the metal af- 
fects electron motions in the planes perpendicular to its 
direction. The effect becomes significant in high mag- 
netic fields when the cyclotron frequency of electrons Q 
exceeds their scattering frequency 1/r. Under such con- 
ditions, so called windows of transparency appear in the 
q, lo space (q, lo are the wave vector and the frequency 
of the electromagnetic wave, respectively). These are re- 
gions where the above described collisionless absorption 



of the external electromagnetic field cannot occur. As 
a result, various kinds of weakly attenuated electromag- 
netic waves whose frequencies and wave vectors belong 
to the regions of transparency (such as helicoidal, cy- 
clotron and magnetohydrodynamic waves) may occur in 
metals 0, All these waves could be excited by an 
external electromagnetic field, and they are generated 
by collective motions of conduction electrons. At low 
temperatures T (Ml 3> kT, k being the Boltzmann's 
constant) the quantization of the electron motions in the 
planes perpendicular to the field B appears. This gives 
rise to oscillations of the electron density of states at 
the Fermi surface (FS) of a metal. These quantum os- 
cillations generate several effects such as de Haas - van 
Alphen oscillations in the magnetization and Shubnikov 
- de Haas oscillations in the magnetoresistivity. Also, the 
quantization of the conduction electrons motions brings 
changes into the geometry of transparency regions, cre- 
ating opportunities for occurrence of particular modes 
which cannot exist otherwise. 

To explain this we assume that a quantizing magnetic 
field is applied along the "z" axis of the chosen coordinate 
system (B = (0,0, B)) . Then the conduction electrons 
acquire Landau energy spectra of the form: 



E(n,p z ,a) = m(n + 1/2) + E {l (p z ,a) 



(1) 



where p z is the electron momentum projection on the 
magnetic field direction, the quantum number n labels 
Landau levels, a is the spin quantum number, and the 
term E\ \ (p z , a) originates from the electron motion along 
the magnetic field. As follows from the Eq. ([1]), at a cer- 
tain value of the magnetic field the "longitudinal" part 
of the energy of an electron at the FS E\ | (p z , a) can 
take on values, belonging to a set of rather narrow in- 
tervals whose widths are of the order of thermal energy 
kT. This leads to a similar distribution of possible values 
of the longitudinal velocity v z = dEi i / dp z at the Fermi 
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surface. Within the low temperature limit (T — » 0), 
these intervals are reduced to the points, and we obtain 
a discrete set of possible v z values for every given B. 
Let us assume that the wave travels along the magnetic 
field (q = (0, 0, q)j, and the condition qv z /Q < 1 is 
satisfied. Under this condition, the charge carriers dis- 
placements in "z" direction during the cyclotron period 
27r/f2 are much smaller than the wavelength A = lit /q. 
Therefore, quasiparticles are seeing nearly the same elec- 
tric field after each turn over their cyclotron orbits. They 
may efficiently absorb the wave energy provided that at 
least one among the permitted values of the longitudinal 
velocity equals the phase velocity of the wave. Due to the 
fact that quantization of the charge carriers motions puts 
essential restrictions on the possible values of v z , extra 
regions of transparency appear thus creating opportuni- 
ties for the specific modes to occur. Possible appearance 
of such weakly damped modes in quasi-isotropic metals 
was predicted by Konstantinov and Perel [H, and they 
were named "quantum waves" [5|. Subsequently, their 
spectra and coupling with ultrasonic waves travelingin a 
metal were theoretically analyzed in several works [6r|8|. 
However, this analysis showed that experimental observa- 
tion of quantum waves in conventional three-dimensional 
metals should be very difficult, because extremely low 
temperatures are required for their appearance. 

In the last three decades various layered conducting 
materials were synthesized. These are organic conduc- 
tors belonging to the family of tctrathiafulvalene salts, 
dichalcogcnidcs of transition metals, graphite and its in- 
tercalates, cuprates and some other. Strong anisotropy 
of the electrical conductivity is inherent for all these ma- 
terials. The anisotropy originates from the fact that con- 
ducting layers are weakly coupled to each other. Corre- 
spondingly, a charge carrier energy is nearly independent 
of the momentum perpendicular to the layers, and the 
charge carriers spectra are quasi- two-dimensional (Q2D). 
Fermi surfaces of Q2D conductors could be described as 
systems of weakly corrugated cylinders 0- 11 1 . These ma- 
terials attract significant interest of the research commu- 
nity, and their electron characteristics are being inten- 
sively studied both theoretically and experimentally. In 
particular, it was shown that the specific FS geometry 
inherent for Q2D conductors provides better opportuni- 
ties for occurrence and observation of weakly attenuated 
electromagnetic waves, which could appear in strong but 
nonquantizing magnetic fields fl2T |. In the present work 
we theoretically analyze spectra of longitudinal "quan- 
tum waves" in Q2D metals and their interaction with 
sound waves. 



II. Dispersions of longitudinal quantum waves 

In our analysis we assume that the field B is orthog- 
onal to the conducting layers, and we use a simple tight- 



binding approximation for the charge carriers dispersion: 



2m± 



2tcos(k z d) (2) 



Here, p± is the momentum projection on the layer plane 
and m± is the effective mass corresponding to the mo- 
tion in this plane, respectively. The parameter t in 
the equation ^ is the interlayer transfer integral, and 
Po = Tth/d where d is the distance between the adja- 
cent layers. In a quantizing magnetic field, the charge 
carriers spectra are described by the Eq. (TTJ) with the 
longitudinal energy of the form: 

Ell {p z ,a) = -2tco S (^A+^ B (3) 

where g is the spin-splitting coefficient and /3o is the 
Bohr's magneton. 

In general, spectra of weakly attenuated waves in met- 
als are determined by the dispersion equation of the form 

GLQ: 



Det 



q 2 Si 3 - q t qj — ^(a^q) 



= 0. 



(4) 



Here, q%,qj are the components of the wave vector q 
within the chosen coordinate system, Sij is the Kro- 
necker symbol and the dielectric permittivity tensor 
£jj (w,q) is related to the conductivity tensor tJij(u), q) : 



^(w,q) = S i:j + — cry(w,q) 



(5) 



So, the last term in the dispersion equation describes the 
response of a metal to an electromagnetic disturbance. 
The electric conductivity tensor is defined by the relation: 



(6) 



where 3 qu> and E qu are the Fourier transforms of the 
charge carriers current density and the electric field, re- 
spectively. To compute the current density, one needs 
to find the expressions for matrix elements of the charge 
carriers density matrix. When the external disturbance 
of the conduction electrons system is weak, the density 
matrix includes the equilibrium term and small nonequi- 
librium correction p, describing the linear response of 
the system to the disturbance. The density matrix p 
satisfies the equation: 



(7) 



where H = H + Hi is the single-particle Hamiltonian 
for charge carriers exposed to an external disturbance. 
In the considered case, the Hamiltonian Hq eigenvalues 
E v are the energies given by the Eqs. (dJ) , (J3J) , and the 



3 



term Hi describes the effects of the electromagnetic dis- 
turbance and scattering. The kinetic equation (J6j> may 
be rewritten in the form 13|, [lj] : 



-{E u , -E u + tvjj)p vv , + - fw)W vvl = I v 
in in 



[p] (8) 



where /„ is the Fermi distribution function for the en- 
ergy E v and I uu i [p] is the collision integral which de- 
scribes the charge carriers scattering. Considering a high 
frequency disturbance (wt > 1) we may use the simple 
approximation: 



Lu> [p] 



-~Pw> 
T 



(9) 



For simplicity, we omit from consideration the effects 
arising due to interactions between the charge carriers. 
Then the term W vv > solely corresponds to the effect of 
the electromagnetic disturbance and can be presented as 
follows: 



Wuu' = en W '(-q)$q W H j vv < 

c 



-q)A 



(10) 



where en vv t{— q) and ] vv i{— q) are the Fourier trans- 
forms of matrix elements of the operators of charge and 
current densities, respectively, and $ qw and A qw are 
the Fourier transforms of the scalar and vector poten- 
tials of the electromagnetic field disturbing the conduc- 
tion electrons. The current density 3 qu equals: 



Jgu = ^Pvv>{u,a L )s u . u (q). 



(11) 



Solving the kinetic equation (jHJ and substituting the re- 
sults for /9 W '(w,q) into Eq. (II ip we may find expressions 
for the conductivity components. 

Below, we consider longitudinal modes traveling along 
the magnetic field so the dispersion equation (|4|) is re- 
duced to the form: 



e zz (io,q) = 0. 



(12) 



Employing Eqs. (|8j)- (|lip and the continuity equation for 
the electric charge: 

g 

— {E v , - £„)n w .(-q) = iqw(-q) (13) 
in 

we may get the following expression for relevant conduc- 
tivity component: 



: ('4)£ 



(U - f V ')n vv '{-q)n v > v {q) 



q A \ lot / *■ — ' E v i — Ey + hcu + ih/r 

vv' 

(14) 

As shown in the Appendix the expression for a zz may 
be converted to the form: 



a zz (u),q) 



=a(q)u{iru + T + + T_ - i(l + ?7 + + ?7_ + A+ + A_)} 

(15) 



where 



*(q) 



irNe 2 m±vo 



CO 

qvo' 



(16) 



A is the mean cross-sectional area of the FS, vq — 
2irt/po is the maximum value of the charge carriers veloc- 
ity component along the magnetic field, ± correspond to 
the different spin orientations with respect to the mag- 
netic field and N is the charge carriers density. The 
Eq. ([15)) includes both real and imaginary parts. It 
gives a good approximation for the considered conduc- 
tivity component at small values of u (it -C 1). The real 
part presents the conductivity as such whereas the imag- 
inary part describes the phase shift of the electric field of 
the incident wave. 

The oscillating terms r\ a and Y a are described by the 
real and imaginary parts of the function X a which has 
the form: 



1 

2qR^ 

r— 1 



Awrt 



(-ir 



D(r) exp 
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X < COS 



Ml 
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U- 



■ COS 



Anrt 

mi 



(17) 



Here U± = (u ± irhq/2po), HHq = Hf3oB is the spin 
splitting energy, R = qv$/Q, F — cA/2irhe. The damp- 
ing factor D{r) = RT{r)Rr>{r) describes deterioration 
of the magnetic oscillations due to the effect of tempera- 
ture {Rx{r)) and charge carriers scattering (Rn{r)). 

Additional oscillating corrections A CT , which appear 
in the Eq. (fT5|) , are closely related to the magnetic quan- 
tum oscillations of the charge carriers density of states 
on the Fermi surface. When the latter is noticeably 
warped (t/Ml > 1) and the cyclotron quantum Ml is 
much smaller than the charge carriers chemical potential 
p (uy p/Ml > 1), the terms A a may be approximated 
as follows 1151: 



A„ = 



(-ir 



2-kH ^ 

r—1 



D(r) 



cos 



2irrF ma 
B 

1 ± min 



D 



z 9 

9 °-0 . 



(18) 



In this expression, i?_o(r) = exp(— nr/flr), so the ef- 
fect of scattering may be accounted for by replacing 
the temperature T in the expression for Rt{t) = 
ry I sinh(r?/) {y = 2ir 2 kT/hQ) by an "effective temper- 
ature" T* =T + T D where T D = h/2irkT is the Dingle 
temperature (l6| . 

Now, we compare two terms in the expression ([5]), and 
we conclude that in the case of e zz , the second term 
significantly exceeds the first one provided that co <C 
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ujq. We omit the lesser term and reduce the dispersion 
equation (|12p to a very simple form: 



Y+ + Y_ = 



(19) 



exp(— 7rrirf / Mlql) . Then we obtain 



(22) 



When the FS becomes perfectly cylindrical (t — > 0), the 
difference between the terms included in the Eq. (fl7|) dis- 
appears, and the oscillating functions rj a and T a turn 
zero. This agrees with the results of a recent work fl5j j . 
where the analysis was carried out in the limit of small q. 
So, the dispersion equation (fT9|) does not have solutions 
corresponding to weakly attenuated waves when 2D con- 
ductors are considered. This conclusion seems reasonable 
because electromagnetic waves in metals are generated 
by collective motions of conduction electrons. At t = 
the charge carriers cannot move between the conducting 
layers of a layered conductor, and this prevents the oc- 
currence of a longitudinal collective mode propagating 
along the perpendicular to the layers. Such modes could 
appear in Q2D metals with moderately corrugated Fermi 
surfaces, where the inequality t > Ml may be satisfied in 
strong magnetic fields. To prove this statement, we ana- 
lyze the behavior of the oscillating terms r\ a and r CT . 

As discussed in the introduction, the quantization of 
the charge carriers motions puts restrictions on the al- 
lowed values of the longitudinal component of the their 
velocity v z = dE/dp z at the Fermi surface. In the low 
temperature and collisionless limit, the v z values at a 
certain magnitude of the field B belong to a discrete set. 
For a Q2D FS determined by the Eq. ([2]), this set consists 
of two subsets related to the FS cross-sections with the 
maximum {u n ^ a } and minimum {u m a } cross-sectional 
areas. Using the Eqs. ([T]), © which determine the charge 
carriers spectra in Q2D conductors in the presence of a 
quantizing magnetic field, we may derive expressions for 
the permitted longitudinal velocities of the quasiparticles 
at the Fermi surface: 



where b = 



t 



t V B 



n 2 



(l-^/l-< CT ), (20) 
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t V 2 



1 F„ 



B 



2 1 



1 - v 2 



(21) 



where w n ^ a = u n , a /vQ and v m ^ — it m ,o-/t>o- It is shown 
in the Appendix that the functions r\ a diverge when U± 
approaches one of the dimensionless velocities introduced 
by Eqs. (|20|) . ([21]) . This result indicates that the disper- 
sion equation p^|) may have solutions corresponding to 
the considered quantum waves. 

Under practical conditions, the combined influence of 
temperature and charge carriers scattering leads to the 
smoothing over the divergencies. To estimate their effect 
on the r\ a oscillations, we again use the well known ex- 
pression Rt{i") — ryj sinh(ry) where y = 2ir 2 kT /Ml. 
Also, we approximate the scattering factor as Rd{t") — 



kT wr 



and £ is Riemann's zeta func- 



tion. As follows from Eq. ([22]), at ujt > u^Jt/kT > 1 
the magnitude of the oscillations is determined by the 
temperature (|?7ct| ~ uMl/kT) whereas at 1 < ujt < 
u^Jt/kT the effect of charge carriers scattering predom- 
inates ( |77cr | ~ qlMl/t). Now, we combine these results 
keeping in mind that oscillating terms given by Eq. (|17p 
must be taken into consideration when the Fermi sur- 
face of a layered conductor is moderately corrugated 
(t > ujpo/q). This follows from the results presented in 
the Appendix. So, we conclude that the function r\ a may 
reach values exceeding unity at the peaks of its quantum 
oscillations provided that: 



Ml 

U fcT >1; 



Ml 
ql- > 1. 



(23) 



where / = vqt. When these conditions are satisfied, the 
dispersion equation (|19l) has solutions corresponding to 
the quantum waves. As follows from the dispersion equa- 
tion, the imaginary terms in the Eq. (|2ip determine the 
quantum waves attenuation. 

The behavior of the functions r\ a and in the low 
temperature and weak scattering region defined by the 
inequalities (|2Tj) is shown in the Fig. 1. In plotting the 
curves presented in this figure we used the relevant data 
reported for a typical Q2D organic metal j3 — {ET^lBr^ 
111 ] . namely: A = 1.31 x lO^^fcgms -1 ) 2 , m± — 
4.5mo (too being a free electron mass) and (v4 max — 
A m in)/A ~ 0.04. Using these results we estimate the 
velocity vq and the transfer integral t, and we get 
vo ~ 4xl0 3 ms _1 and t ~ 2meV, respectively. As shown 
in the Fig. 1, contributions to r/ a from both maximum 
and minimum cross-sectional areas of the Q2D Fermi sur- 
face reach values exceeding unity at the peaks of oscil- 
lations at T~ lOmK. Presently, such temperatures are 
quite accessible in experiments. As for the term r CT , its 
magnetic field dependence reveals a sequence of nearly 
rectangular peaks slightly smoothed near the vertexes 
due to the effects of temperature and scattering. The 
peaks heights and widths are determined by the value 
of the parameter qR. In the limit of long wavelengths 
(qR <C 1), the peaks become very high and narrow sim- 
ilar to those describing so called giant quantum oscil- 
lations in the ultrasound attenuation observed in some 
conventional 3D metals [l|. In the present case, these 
peaks appear as a result of collisionless absorption of the 
electromagnetic waves energy by charge carriers. The 
effect of giant quantum oscillations of electromagnetic 
field attenuation was predicted for quasi-isotropic met- 
als in the works 17], |l8j. The peaks are separated by 
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FIG. 1: (Color online) Magnetic quantum oscillations in the 
real (left panel) and imaginary (right panel) parts of the func- 
tion Y a . Solid and dashed lines correspond to the contribu- 
tions from minimum and maximum cross-sections of a mod- 
erately corrugated Fermi surface of a Q2D conductor. The 
curves are plotted assuming that F a /Bo = 100, t/M}(Bo) = 
4, T = WmK, B = 20T, m± = 4.5m . 



FIG. 2: (Color online) The dependencies of the real (left 
panel) and imaginary (right panel) parts on a zz of uj/qvo. 
The curves are plotted neglecting the spin-splitting of the 
charge carriers energy levels to simplify the shape of the peaks 
in <r zz real part. However, the contributions form both max- 
imum and minimum cross-sectional areas are included. The 
curves are plotted assuming B = Bo — 20T and using the 
same values for the remaining parameters as in the Fig. 1. 



the intervals where T a equals —ttu/2 and appears com- 
pletely balanced by the semi-classical contribution to the 
imaginary part of Y a . 

To further clarify the physical meaning of Y a we turn 
to the expression (fl"4|) for a zz (uj,q). We remark that 
the condition of vanishing of the imaginary part agrees 
with the dispersion equation for the undamped quantum 
waves. The dependencies of the real and imaginary parts 
of a zz on the frequency w at a certain fixed value of 
the magnetic field are shown in the Fig. 2. The real 
part exhibits peaks within some frequency intervals and 
remains zero beyond them. Within these regions of high 
conductivity charge carriers intensively absorb the energy 
of the incident wave preventing its propagation in the 
metal. The peaks in the conductivity appear as a result 
of summing up contributions from T + and F_. At low 
temperatures both T a reveal sequences of nearly rect- 
angular peaks. The relative positions of nearest peaks 
in r + and T_ are determined by the magnitude of the 
magnetic field B, the value of the transfer integral t, 
which controls the difference between the maximum and 
minimum cross-sectional areas of the FS, frequency u> 
and spin splitting characterized by g^f ■ Varying these 
parameters one may observe subsequent variations in the 
relative positions of the nearest peaks. At some values 
of the relevant parameters th nearest peak T+ and F_ 
are completely or partially overlapped, at order values 
they stay apart from each other. As a result, the con- 
ductivity dependence of frequency may reveal moderately 
low trapezoid-like peaks corresponding to separate near- 
est peaks in T a for different spin orientations or higher 
features of a characteristic bottle-like shape, which ap- 
pear due to the partial overlapping of the nearest peaks 
in r a . Such features are shown in the Fig. 2. 



High trapezoid-like peaks occurring, which the nearest 
peaks in r + and T_ completely overlap may appear in 
the frequency dependence of Re a zz , as well. Longitudi- 
nal quantum waves may propagate in the regions of low 
conductivity where the metal becomes "transparent" for 
the electromagnetic fields. The locations of the regions 
of transparency are determined by relations: 

2 , t 7" 2 - 2 2 ^- j r2 , 2 

W <U + < w na ; w n , +h(7 , <U_ < w n ,^,\ 

(24) 

These inequalities may be consistently satisfied when 
Landau quantum numbers n,n' ' ,m,m! and spin num- 
bers a, a 1 accept fitting values. They specify the regions 
where the dispersion equation (fl~9|) has solutions corre- 
sponding to weakly attenuated longitudinal modes prop- 
agating along the magnetic field B. 

19 Now, we analyze the effect of the corrections on 
the dispersions of the quantum waves. The magnitude of 
A CT at the peaks of oscillations may be estimated as: 

f mi nn\ 1/2 

In principle, A CT magnitude may reach values compared 
to unity when the effective temperature T* becomes ex- 
tremely low thus providing negligibility of thermal and 
scattering effects. However, this may happen at signif- 
icantly lower temperatures and weaker scattering than 
those needed for the terms r) a to reach values of the 
order of unity at the oscillations peaks. To justify this 
statement, we take the same values for the parameters 
characterizing the Fermi surface which were used in plot- 
ting of the Fig. 1, namely: B = 20T, u 10 12 s -1 , 
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ujt ~ 1, and T = lOmK. This gives for the Dingle 
temperature the estimation ~ lOOmif, thus provid- 
ing that peak values of A CT are approximately ten times 
smaller that the corresponding values of r\ a and signif- 
icantly smaller than unity. On these grounds, we omit 
the corrections A CT in further analysis. 

20 Nevertheless, it is worthwhile to remark that the 
charge carriers scattering differently affects r\ a and A CT . 
This difference reflects the diversity in the nature and ori- 
gin of these oscillating contributions to Y a . The terms 
A CT appear solely due to quantization of the charge car- 
riers motion in the planes perpendicular to the strong 
magnetic field. The charge carriers motions along the 
magnetic field direction do not noticeably modify these 
terms. Naturally, the effect of scattering is completely de- 
termined by the value of fir, which indicates how long a 
charge carrier could stay at a certain cyclotron orbit. On 
the contrary, charge carriers motions along the magnetic 
field are crucial for occurrence of the terms 7] a . There- 
fore, the effect of scattering is mostly specified by the 
value of ql which indicates how long a charge carrier 
could preserve its motion "in phase" with the electro- 
magnetic wave absorbing its energy. 

21 We can easily solve the dispersion equation (fT9|) 
within the long wavelength limit qR <c 1. In this case, 
the inequalities ([24| determining the regions of trans- 
parency for the considered waves, may be reduced to the 
form: 



< u < m 



V, h+ < U < V n+ i.- 



(26) 



where ± indicate two spin orientations. The terms r\ a 
given by Eq. (I17[) may be approximated as follows: 



T] a = — U < cot 



2tt 



cot 



2tt 



max erg Oq t 

b 2 n Ki 



(27) 



In the expression (|27[) . the cotangents diverge when u 
approaches one of the velocities ai„ jff and/or v nM . As 
follows from Eqs. (fTTJ)) , (|2"7|) . within the long wavelength 
limit, the weakly attenuated longitudinal modes deter- 
mined by the Eq. (|19|) . have linear dispersions, and their 
speeds strongly depend of the magnetic field magnitude, 
being close to w n>a or v n _ a . After some simple algebra, 
the dispersion equation for the mode occurring in the 
window of transparency specified by the first inequality 
(f26|) . accept the form: 



1 



W n +1,- 



4f wi 



4i wl_ 



= 



(28) 



Similarly, we obtain the dispersion equation for the mode, 
which may appear in the window, determined by the sec- 
ond inequality (|26|) : 



Ml 



At u 



At u 2 



0. 



(29) 




Ql LJ , , 1 Ql , , 1 1 1 

ff.95 1 1.05 1.1 1.15 1.5 1.55 1.6 1.65 1.7 
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(B n /B-i)xi(r 



FIG. 3: (Color online) Magnetic field dependencies of slow 
longitudinal quantum modes propagating along the magnetic 
field related to the maximum (left panel) and minimum (right 
panel) cross-sections of a Q2D Fermi surface. The speeds of 
waves shown in the figure are close to the reduced longitudinal 
velocities of the charge carriers w n ,+ , (left panel) 

and v n ,+ , Vn+i,- (right panel). The curves are plotted using 
the same values of relevant parameters as in the Fig. 1. 



The solutions of these equations corresponding to the 
weakly damped longitudinal modes are shown in the Fig. 
3. The curves plotted in this figure represent modes re- 
lated to the maximum and minimum cross-sections of the 
Fermi surface, respectively. They are picked out of the 
sets of branches, each branch corresponding to a certain 
value of the Landau number. 

Low frequency modes, whose speeds have the same or- 
der as the speed of sound waves propagating along the 
axis of the Q2D Fermi surface may be excited by the 
latter. In the next Section we consider the interaction of 
quantum waves with the sound and its effect on quantum 
oscillations of the sound speed. 



III. Interaction of the longitudinal quantum 
waves with sound waves 



Elastic response of a metal to an external deforma- 
tion includes rearrangements in the system of conduc- 
tion electrons. When a sound wave travels in the metal, 
it gives rise to a collective motion of the charge carriers 
producing a self-consistent electromagnetic field accom- 
panying the wave. Also, the lattice deformation gives 
rise to corrections to the crystalline fields, which could be 
described by means of deformation potentials. However, 
in the considered case accounting for these corrections 
does not bring any qualitative changes in the results, so 
we omit them for brevity. The electron contribution to 
the elastic response of the metal manifests itself in mag- 
netic oscillations of sound velocity and attenuation, and 
in the coupling of electromagnetic waves, occurring in the 
metal in the presence of an external magnetic field, to the 
sound waves of fitting polarizations and propagation di- 
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rcctions. Here, we consider interaction of the longitudinal 
quantum waves traveling along the symmetry axis of the 
Fermi surface of a Q2D conductor with the longitudinal 
sound waves propagating in the same direction. Apply- 
ing general equations for the magnetoacoustic response 
of a metal LJ, [l4 1 , we may write the following equation 
for the lattice displacement vector u; = (0,0, u/) : 



Pr, 



d 2 ui d 2 ui 



dt 2 



= X- 



dz 2 



QE 



(30) 



where p m , A are the material density and the appropri- 
ate elastic constant of the lattice, respectively; Q is the 
charge density of the lattice, and E = (0, 0, E) is the 
electric field accompanying the lattice displacements due 
to the charge carriers response. The field E obeys the 
Maxwell equations, which gives: 



1 d 2 E _ 4tt 

— *J Z 



dt 2 



(31) 



Here, the term Qdui/dt corresponds to the lattice con- 
tribution to the electric current, whereas J z (a,t) is the 
contribution from the charge carriers. Now, we assume 
that the time and space dependencies of the lattice dis- 
placement ui and the electric field E have a simple 



harmonic form, namely: ui 



, exp(iqz — iuit), E 



E qul exp(iqz — iu)t). As follows from Eqs. (|30)) . (|3T 



E. 



quj 



q 2 Qu qiu 

e 2 9o 



(32) 



where go is the charge carriers density of states on the 
Fermi surface in the absence of the magnetic field and Y± 
are introduced by the Eqs. (15]). Using Eqs. ([30"]) -(f32" l) 
we may derive the dispersion equation for the coupled 
quantum waves and the sound: 



2 2 , 

q s + 



q 2 N 2 



p m g {Y+ + y_) 



(33) 



In this equation, sq = \/A/p m is the "lattice contribu- 
tion" to the speed of sound, and N is the charge carriers 
density. The equation (|33|) agrees with the correspond- 
ing result obtained for the particular case of an isotropic 
metal with a spherical Fermi surface [8| . From this equa- 
tion one could see that when the denominator of the 
second term on the right hand side turns zero indicat- 
ing the quantum wave appearance, this strongly affects 
the sound frequency. Thus, the dispersion equation (|33|) 
gives a positive evidence of interaction between the sound 
and quantum waves. 

We can analytically solve Eq. (|33|) in the long wave- 
length limit qR <C 1 provided that the spin splitting 
energy ^MIq equals the cyclotron quantum Ml. In this 

particular case the longitudinal velocities of the charge 
carriers associated with two subsequent values of the Lan- 
dau number and different spin orientations coincide. For 



certainty, we assume u>n,+ = w n +i,- = w n , and we an- 
alyze the solutions of the Eq. (|33| with u close the w n . 
We approximate: 



7TU 

V±~ — cot 



nt I 2 



u 2 ) 



At w 2 — u 



(34) 



Using this approximation, (|33[) could be reduced to a 
biquadratic equation, whose solutions have the form: 



1 



1 



Here, «o — s/v , s = ysg + s\ is the speed of sound 
including the contribution from the charge carriers s\ = 
\J N 2 /p m ga, and U n is the speed of the quantum wave 
determined by the expression: 



^ n 



w„ 



Ml 
~2t 



-Wr, 



(36) 



The dispersion equation solutions of the form (|35j) are 
shown in the left panel of the Fig. 4. The plotted curves 
obviously correspond to the coupled sound and electro- 
magnetic quantum waves of the same polarization and 
propagation direction. In more practical situation when 
^Mlo 7^ hil, the dispersion equation (|33[) describes cou- 
pled quantum and sound waves, as well. 

The interaction with the quantum waves makes a sig- 
nificant effect on the magnetic quantum oscillations of the 
the sound velocity. When the temperature is moderately 
low (u~ x > Ml/kT > 1) and the charge carriers scatter- 
ing is moderately weak (1 < ql < t/Ml), the oscillations 
of the functions r\ a are small in magnitude. As shown in 
Ref. [H|], in this regime the oscillating correction to the 
velocity of sound includes two terms, one being propor- 
tional to r) + + ?7_ and another one to A + + A_. These 
terms describe magnetic oscillations different in shape, 
phase and period but their magnitudes are of the same 
order. However, at very low temperatures and weak scat- 
tering, when the inequalities (|23p are satisfied, the peak 
values of the functions r\ a may reach and exceed unity. 
At first glance, this gives grounds to conjecture that in 
this case the magnetic quantum oscillations of the sound 
velocity would be significantly increased in magnitude. 
However, in the regime specified by Eqs. (|23j) . the propa- 
gation of the sound wave is accompanied by generation of 
the sequence of quantum waves branches. This is shown 
in the right panel of the Fig. 4. Tracing the sound branch 
in this figure, we see that jumps in the sound velocity 
related to the interaction with the quantum waves are 
rather small. Their magnitudes significantly concede the 
sound velocity value s taken in the absence of an exter- 
nal magnetic field. 

Finally, we briefly consider magnetic field dependence 
of the longitudinal component of the charge carriers spin 
density s z = s quj exp(iqz — iojt) under the conditions 
determined by Eqs. J23[) . The spin density amplitude 
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IV. Concluding remarks 
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FIG. 4: (Color online) Coupling of slow longitudinal quan- 
tum waves to the longitudinal sound propagating along the 
magnetic field in a Q2D metal. The quantum mode is related 
to the maximum cross-sectional area of the Fermi surface. 
The curves are plotted assuming that the relevant parame- 
ters accept the same values as in the previous figures. 



Squ is proportional to the difference between the densities 
of the charge carriers with different spin orientations, so 
we may write: 



-hr,+ - ~\ 



(37) 



Here, the Fourier transform of the charge carriers density 
is described by the expression: 



n qu = _5_3/w(w,q)jVi,(q) 

77' 



(38) 



where 7 is the set of the orbital quantum 
numbers (-{V} = {7,cr}) and the density matrix p 
is detrmined by the kinetic equation (|SJ). As follows 
from (|38|) . the charge carriers density Fourier transform 
is proportional to E qul : 



n qul = E qbJ Y±. 



(39) 



Combining this result with the expression (|32l) we get: 



nf u = inqQu qbJ 



y+ + y_ 



So, we obtain: 



— iqU ' 

2 Y+ + Y- 



(40) 



(41) 



This shows that the quantum waves generated by com- 
bined effects of the lattice vibrations and the collective 
motions of charge carriers providing the electroneutral- 
ity of the whole system, are accompanied by magnetic 
oscillations of the longitudinal spin density. 



In the present work we theoretically analyzed longitu- 
dinal collective modes which could propagate along the 
quantizing magnetic field in Q2D conductors. Within the 
chosen geometry the magnetic field was directed at the 
right angle to the conducting layers. We showed that 
these "quantum waves" may occur at low temperatures 
and weak charge carriers scattering provided that the 
Fermi surface warping is rather pronounced (t > ujpa / q) . 
Especially interesting is the case of slow quantum waves 
whose speeds are close to the speed of sound propagating 
in a Q2D metal perpendicularly to the conducting planes. 
Such slow modes could interact with the sound waves, 
and this interaction significantly affects magnetic quan- 
tum oscillations of the sound velocity. As was mentioned 
above, the possible appearance of longitudinal quantum 
waves in conventional quasi-isotropic metals was theo- 
retically predicted a few decades ago but they were not 
observed in experiments. Q2D conductors provide much 
better opportunities for such observations. Now, we 
briefly discuss which properties of Q2D conductors give 
them the advantage over conventional metals in this re- 
spect. The conditions on the temperature and the inten- 
sity of charge carriers scatte ring in quasi-isotropic met- 



als were presented in Ref. 14(. They could be writ- 
ten in the form (|23p provided that transfer integral t in 
the second inequality is replaced by the charge carriers 
chemical potential fi. The restrictments put upon the 
scattering intensity does not significantly diverse for 3D 
and Q2D materials for the parameter po usually takes 
on values of the same order as the Fermi momentum in a 
quasi-isotropic metal. However, a pronounced difference 
reveals itself in the temperature regime required to pro- 
vide the quantum waves occurrence in Q2D and 3D met- 
als. Due to the specifics of the charge carriers dispersions 
in Q2D materials the maximum charge carriers velocity 
in the direction perpendicular to the conducting layers is 
much smaller than typical Fermi velocities of conduction 
electrons in quasi-isotropic metals. Assuming that quan- 
tum waves are traveling with speeds close to the speed of 
sound propagating along the same direction, the param- 
eter u in Q2D conductors probably takes on values no 
less than 0.1, whereas in 3D metals it has the order of 
10~ 3 — 10~ 4 . Therefore, within the chosen geometry the 
considered quantum waves may occur in Q2D conductors 
at temperatures significantly exceeding those required for 
them to appear in 3D metals. As shown in this work, one 
may expect to observe the longitudinal quantum waves 
traveling along the magnetic field in Q2D organic metals 
and their effect on the magnetic quantum oscillations of 
the sound velocity at quite accessible temperatures of the 
order of lOmK. Such observations should provide better 
insight in the specific features of low temperature mag- 
netic quantum oscillations of the elastic response of Q2D 
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conductors. 
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V. Appendix 

Expressions for the functions Y± , which appear in the 
Eqs. (fT5")) and (fT9l are given by: 




p Rez 



1 



1 



: 47r 2 SA 2 5o E 



dp z 



f(n,p z ,a) - f{n,p z - hq,a) 

E(n,p z — hq, a) — E(n,p z ,<r) + frw + ih/i 



(42) 



where A is the magnetic length. 

Using the Poisson summation formula 



E 

71 = 



ip [n + 



1 



r— — oo 



(f(x) exp(27rzrx) 



(43) 



we may present the functions Y a as sums of terms Yo a = 

l^t^o — (l + = ~(l + independent of the 

2TT 2 h 3 gp\ 2 J 2 V 2 / 

magnetic field magnitude and oscillating contributions 
Y a . The latter could be approximated by expressions: 



Y a = Y a (q)+Y a {-q) 



(44) 



where 



-I OO r. OO 

( 2 ™r g Jo 

cos[2ttt F(E,p z , a) I B]dp z 



_ PQ E(n,p z ± Tig, a) - E(n,p z ,a) T/wTft/r 



(45) 



FIG. 5: (Color online) The schematic of the contour chosen 
to carry out integration over p z in the expression (|45[) . The 
point indicates a pole in the integrand. 



where u — uj/qvo, I — vqt, and vo — 2nt/po is the 
maximum value of the charge carriers velocity component 
along the magnetic field. The poles are situated within 
the contours of integration when the Fermi surface corru- 
gation is sufficiently pronounced to satisfy the inequality 
t > ujpo/q. The residues from these poles give contri- 
butions to the oscillating parts of the functions Y„ . Be- 
sides, the oscillating parts include corrections originating 
from the integration over the semicircular arcs. Assum- 
ing that u 1 we obtain the following approximations 
for Y a : 



VKU 



iT r 



(48) 



Here, r\ a and T s igma are equal to the real and imag- 
inary parts of the function X a , respectively. The func- 
tion X a has the form (|17p Neglecting for a while the 
effects of temperature and scattering, we may carry out 
summation over "r" in the expression (|17p . using the 
formula 11911 : 



( -l)^-^ln{2cos-}. 

fc=i 



(49) 



Here, 



F(E, Pz ,a) 



2-nhe 



A(E, Pz ,a) 



(46) 



and A(E,p z ,a) is the cross-sectional area of the corre- 
sponding isoenergetic surface. 

To properly estimate the integrals over p z in Eq. (|45|) 
we extend integrands over the upper halves of the com- 
plex plane and we choose the integration paths consisting 
of segments of the real axes — po < p z < Po and semicir- 
cles of the radii po (see Fig. 5). The integrands possess 
poles, whose positions are determined by equations: 



Pa 



irhq I p 

± COS 7T 

2po 



Po 



u- — =0 
ql 



(47) 



We get 



cos 



^ = 4^? ln 



+ In 

AqR 



cos 



cos 



■Frr 



B 



(jg^lo 





u 2 



1 



F„ 



B 



B 



n 

agQo 



l 



B 



(50) 



where F m i n and -F max correspond to the minimum and 
maximum cross-sectional areas of the Fermi surface. 

When U± approaches one of the dimensionless ve- 
locities given by the Eqs. (|2"0|) . (f2~Tj). the value of the 
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corresponding cosine in the Eq. (|50|) approaches zero, 
and the function r\ G diverges. For weakly corrugated 
FS (t < ujpo/q), the pole is situated without the con- 
tour of integration, therefore the terms r\ a and T s igma 
disappear from the expression for Y s igma. Besides the 
contributions from the pole, the expression (|48p includes 
the terms A CT originating from the integration over the 
semicircular are shown in the Fig. 5. At t/fkv > 1 these 
terms may be presented as follows 14|: 



i / m ^i-iy 



(r){ 



exp 



2nirF£ lax 



B 



T 



G 



exp 



B 



ITT 

T 



(51) 



where 



G 



i 

± 2 



exp[±iy] exp 



2nJrF° 



i . u 2 v/B dy. 

(52) 

The ratio F/B has the same order as fi/KD, where \x 
is the charge carriers chemical potential and h£l is the 
cyclotron quantum. Under the condition u\J fi/hfl > 1 
we may approximate the functions using asymptotic 
expressions for Fresnel integrals, and we get Gf ps 1. 
Substituting these approximations into the Eq. (I5ip we 
obtain: 



r—1 v 



2?rrF„ 



D(r) 



cos 



B 

ItTT 'Frnin 

B 



9 

irra — — 
2 fi 

"TO — + 

2 SI 



(53) 
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